mm m 

MICROBIOLOGY 



ORIGINAL RESEARCH ARTICLE 

published: 20 December 2013 
doi: 10.3389/fmicb.2013.00399 




Patterns of microbial diversity along a salinity gradient in 
the Guerrero Negro solar saltern, Baja CA Sur, Mexico 

Jesse G. Dillon * Mark Carlin, Abraham Gutierrez, Vivian Nguyen and Nathan McLain 

Department of Biological Sciences, California State University, Long Beach, CA, USA 



Edited by: 

H. Thane Papke, University of 
Connecticut, USA 

Reviewed by: 

B. Thane Papke, University of 

Connecticut, USA 

Lejia Pasic, University of Ljubljana, 

Slovenia 

*Conespondence: 

Jesse G. Dillon, Department of 
Biological Sciences, California State 
University 1250 Bellflower Blvd., 
Long Beach, CA 90840, USA 
e-mail: jesse. dillon@csulb. edu 



The goal of this study was to use environmental sequencing of 16S rRNA and bop genes 
to compare the diversity of planktonic bacteria and archaea across ponds with increasing 
salinity in the Exportadora de Sal (ESSA) evaporative saltern in Guerrero Negro, Baja CA S., 
Mexico. We hypothesized that diverse communities of heterotrophic bacteria and archaea 
would be found in the ESSA ponds, but that bacterial diversity would decrease relative 
to archaea at the highest salinities. Archaeal 16S rRNA diversity was higher in Ponds 
11 and 12 (370 and 380 g r'' total salts, respectively) compared to Pond 9 (180g r'' 
total salts). Both Pond 11 and 12 communities had high representation (47 and 45% 
of clones, respectively) by Haloquadratum walsbyi-Wke (99% similarity) lineages. The 
archaeal community in Pond 9 was dominated (79%) by a single uncultured phylotype 
with 99% similarity to sequences recovered from the Sfax saltern in Tunisia. This pattern 
was mirrored in bop gene diversity with greater numbers of highly supported phylotypes 
including many Haloquadratum-Wke sequences from the two highest salinity ponds. In 
Pond 9, most bop sequences, were not closely related to sequences in databases. 
Bacterial 16S rRNA diversity was higher than archaeal in both Pond 9 and Pond 12 samples, 
but not Pond 11, where a non-Salinibacter lineage within the Bacteroidetes >98% similar 
to environmental clones recovered from Lake Tuz in Turkey and a saltern in Chula Vista, CA 
was most abundant (69% of community). This OTU was also the most abundant in Pond 
12, but only represented 14% of clones in the more diverse pond. The most abundant OTU 
in Pond 9 (33% of community) was 99% similar to an uncultured gammaproteobacterial 
clone from the Salton Sea. Results suggest that the communities of saltern bacteria and 
archaea vary even in ponds with similar salinity and further investigation into the ecology 
of diverse, uncultured halophile communities is warranted. 
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INTRODUCTION 

Some of the best examples of chemical gradients are found in 
solar salterns around the world (Anton et al., 2000; Litchfield 
et al, 2001; Baati et al, 2008; Oren, 2008; Manikandan et al., 2009; 
Oh et al, 2010). Like most hypersaline habitats, these evaporation 
ponds typically contain abundant microbial populations includ- 
ing members of all three domains of life (Javor, 1989). Traditional 
cultivation studies as well as molecular sequencing and FISH 
studies in salterns have revealed diverse communities dominated 
by phototrophs like Dunaliella as well as aerobic heterotrophic 
prokaryotes (Anton et al, 1999; Benlloch et al, 2001; Pasic et al., 
2005; Maturrano et al, 2006; Papke et al, 2007; Baati et al, 2008). 
However, biogeographic differences in the specific communities 
found in salterns have been observed, perhaps due to dispersal 
limitation (Oren and Rodriguez- Valera, 2001; Zhaxybayeva et al., 
2013). 

Variation in prokaryotic communities along salinity gradi- 
ents has also been reported as studies have found that microbial 
species richness decreases with increasing salinity, often result- 
ing in a few dominant phylotypes found in highest salinity ponds 
(Casamayor et al, 2000; Benlloch et al, 2002; Baati et al, 2008). 



In many salterns the highest salinities are dominated by haloar- 
chaea such as Haloquadratum and Halorubrum (Anton et al., 
1999; Burns et al, 2004; Maturrano et al, 2006; Oh et al, 2010) 
as well as having significant proportions (up to 15-27%) of 
extremely halophilic members of Salinibacter (Anton et al., 2000; 
Oren, 2002; Ovreas et al, 2003; Anton et al., 2008). However, in 
other salterns, Haloquadratum (Pasic et al., 2005) or Salinibacter 
(Maturrano et al, 2006) are absent or rare. 

The saltern Exportadora de Sal (ESSA) in Guerrero Negro, Baja 
CA S., Mexico covers an area over 300 km^ and is the world's 
largest producer of evaporative salt. Salt water is pumped year- 
round from the modestly hypersaline Ojo de Liebre lagoon and 
than slowly pumped through a series of large (many >lkm^), 
shallow (~1 m deep), interconnected ponds. The ponds display a 
gradient of chemical make-up as salts precipitate (e.g., gypsum) 
as the water moves up the evaporation scale for over a year (Javor, 
1983). As with many salterns, environmental conditions are quite 
stable over time. Nutrient levels are relatively low at the ESSA 
saltern and it has been classified as oligotrophic (Javor, 1983, 
1989). The ESSA saltern has been intensively studied. However, 
most of the research on the diversity of gene sequences of 
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microbes in this system has focused either on the well-developed, 
benthic microbial mats found at moderate salinities (~70-100g 
r^) (Nubel et al, 2001; Ley et al, 2006; Feazel et al, 2008; 
Kunin et al., 2008; Dillon et al, 2009b; Robertson et al., 2009) or 
within evaporites (Sahl et al., 2008). Few studies have been per- 
formed in the ESSA ponds at higher salinities (> 150-160 g 1^'), 
where planktonic communities dominate (Javor and Castenholz, 
1981) and these were primarily cultivation-based (Javor, 1984; 
Sabet et al., 2009). This study represents the first use of culture- 
independent, molecular techniques to examine the diversity of 
planktonic microbes in the ESSA saltern at the highest salin- 
ities where benthic microbial mats are not found. Based on 
the prior cultivation work, we hypothesized that diverse com- 
munities of both halophilic bacteria and archaea would be 
found in the ESSA ponds, but that bacterial community diver- 
sity would decrease in comparison with archaea at the highest 
salinities. 

METHODS AND MATERIALS 

WATER SAMPLE COLLECTION AND ANALYSIS 

In conjunction with a broader study aimed at cultivation of 
halophUes (Sabet et al., 2009), water samples were collected in 
February, 2006, from three ponds along a salinity gradient at 
the ESSA saltworks, Guerrero Negro, Baja California Sur, Mexico. 
Replicate 50 ml water samples (n = 3-5) were collected via near- 
shore surface grabs in evaporative ponds (Ponds 9 and 11) and 
in a crystallizer pond (Pond 12). Samples were frozen in liq- 
uid nitrogen for transport back to CSU Long Beach and stored 
at — 80°C prior to analysis. Pond 9 bottom was covered with gyp- 
sum precipitate, while Pond 12 contained halite with Pond 11 
showing a soft sediment bottom with evidence of both gypsum 
and halite precipitates. Salinity (total salts) of each water source 
was measured using a refractometer and when salinities were off 
scale (i.e., >280g 1^^) by dilution prior to reading. Water tem- 
perature was measured for each pond using a handheld probe 
(Russell RL060P, Thermo Electron Corp., Beverly MA, USA). 
Water samples (n = 1) were filtered through a 0.2 |im filter, 
diluted, and analyzed for major cation content using inductively- 
coupled plasma mass spectrometry (ICP-MS) at the Institute 
for Integrated Research in Materials, Environments, and Society 
(IIRMES) on the CSULB campus and anions via ion chro- 
matography (882 Compact IC plus, Metrohm, Riverview, FL, 
USA) using the EPA 300.0 method at Physis Environmental 
Laboratories (Anaheim, CA). 

DNA EXTRACTION 

The 50 ml water samples were thawed and pre-filtered through 
a lOjxm dia. pore size nylon membrane filters (GE Osmonics, 
Minnetonka, MN) to remove large particles and algae. The water 
was re-filtered through a 0.22 [Lm polysulfone membrane (GE 
Osmonics) to collect the bacterial and archaeal cell fraction. 
Nucleic acids were extracted from the filter using a modified 
protocol of Benlloch et al. (2001). Filters were cut into pieces, 
washed with 2 ml of sterile nanopure water, vortexed, and the 
supernatant treated with SDS (1% w/v) and proteinase K (0.5 mg 
ml^') and the samples incubated at 55° C for 2h, then boiled 
for 2 min. Nucleic acids were extracted twice with 1 volume of 



phenol/chloroform/isoamyl alcohol (lAA) (50:49:1), centrifuged 
at 3300 X g for 20 min, and extracted again with an equal 
volume of chloroform:IAA (49:1) and centrifuged for an addi- 
tional 5 min. The aqueous supernatants were precipitated in 
2 volumes of 100% ethanol and centrifuged at 3220 x g for 
20 min at 4°C. The supernatant was decanted and the tube 
was allowed to dry. The pellet was resuspended in 100 (A L ster- 
ile nanopure water, incubated at 55° C for 30 min and stored 
at -20°C. 

PCR, CLONING, AND SEQUENCING 

PGR amplifications of 16S rRNA were performed using puri- 
fied nucleic acid and bacterial 16S rRNA primers GM3f (5'- 
AGAGTTTGATCMTGGC) and GM4r (5'-TACCTTGTTACGAC 
TT) (Muyzer et al., 1995) and archaeal 16S rRNA primers 
arch21f (TTCGGGTTGATCCYGCCGGA) (Belong, 1992) with 
either the archaea-specific 958r (YCCGGCGTTGAMTCCAATT) 
(Delong, 1992) or the universal 1392r reverse primer (Stahl et al., 
1988). Amplifications of bacteriorhodopsin (bop) genes used the 
bop401F (GACTGGTTGTTYACVACGCC) and bop795R (AAG 
CCGAAGCCGAYCTTBGC) primers (Papke et al, 2007). The 
20 |iL reaction mixtures contained 1 x PCR buffer (Invitrogen, 
Carlsbad, CA), 0.2 mM each dNTP (Promega, Madison, WI), 
10 pmol each primer (Operon, Huntsville, AL), 1 U Platinum Taq 
polymerase (Invitrogen), and 50-100 ng of purified nucleic acids. 
For most samples, 1 \i\ of bovine serum albumin (0.4% w/v) was 
added to reaction mixtures to facilitate amplification. Reaction 
conditions were as follows: initial denaturation (94° C for 5 min) 
followed by 30 cycles of denaturation (94°C for 30 s), annealing 
(53°C for 30 s), and extension (72°C for 90 s) and a final extension 
(72°C for 10 min) using a mastercycler (Eppendorf, Hauppauge, 
NY). The resulting amplicons were ligated into pCR4 TOPO 
vector with the TOPO® Cloning PCR Cloning Kit (Invitrogen, 
Carlsbad, CA) and transformed into One Shot® TOP 10 chem- 
ically competent E. coli cells according to the manufacturer's 
instructions. Transformants were plated on LB plates with 100 |ig 
L^ ^ ampicillin. For each sample, colonies were picked with a ster- 
ile toothpick and grown up in 75 [iL of LB -|- ampicillin broth in 
96-well plates, diluted to a final concentration of 15% (w/v) with 
sterile glycerol, and stored at — 80°C. For all clones, cells were 
grown in 2 ml of LB -|- ampicillin broth at 37°C overnight and 
plasmid minipreps performed using the GenCatch plasmid DNA 
purification kit (Epoch Biolabs, Sugar Land, TX). Plasmids were 
sequenced with M13forward and reverse primers by a commercial 
sequencing facility (University of WA High Throughput Genomic 
Center, Seattle, WA). 

PHYLOGENETIC AND STATISTICAL ANALYSES 

For 16S rRNA a total of 296 archaeal sequences and 254 bacterial 
sequences were obtained from 1-2 clone libraries for each pond 
separately for each gene. Chimera detection of sequences was per- 
formed using Mallard software analysis (Ashelford et al., 2006) 
and Pintail (Ashelford et al, 2005). Non-chimeric 16S rRNA 
sequences were aligned using the SINA aligner on the SILVA web- 
site (Pruesse et al, 2007) and imported into ARB software and 
manually refined with reference to nearest neighbor taxa from the 
V. 102 database (Ludwig et al, 2004). A total of 167 high quality 
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bop sequences were derived from a total of 288 clones generated 
via a single 96-clone library from each pond. These sequences 
were initially aligned to all available gene sequences in the NCBI 
website using Clustal X (Larkin et al., 2007), and then imported 
into a custom-created ARB database. Custom lane masks of 
aligned sequences were created excluding hypervariable regions 
(16S rRNA genes) and ambiguous nucleotide positions. This 
resulted in the export of 874 nt (archaeal 16S rRNA) and 1033 nt 
(bacterial 16S rRNA) and 312 nt (bop). Maximum Likelihood 
trees were constructed using the Blackbox RaxML tool on CIPRES 
Science Gateway v. 7.2 (Miller et al, 2010) with 750 boot- 
strap pseudoreplications or fewer if stopped using the automated 
MRE bootstopping criterion (e.g., 250, bacterial 16S rRNA). 
Where full-length sequences were not obtained, partial sequences 
were added to the 16S rRNA trees using the parsimony tool 
in ARB. 

Additional statistical analyses of the three ponds' clone library 
sequences were performed using distance matrices generated in 
ARB (Ludwig et al., 2004). These datasets differed from those 
used for phylogenetic analyses in that redundant sequences were 
included, and new custom filters were made for archaeal (304 nt), 
bacterial (412 nt), and bop (356 nt) sequences. Exported sim- 
ilarity matrices were used to cluster sequences into OTU by 
pair-wise sequence identity with the average neighbor algorithm 
at a evolutionary distance of 0.00 (actually represents <0.005), 
0.01, 0.03, and 0.05 (Schloss et al, 2009) corresponding to the 
100, 99, 97, and 95% similarity cut-offs, respectively. Rarefaction 
curves were generated using a resampling without replacement 
approach. Estimations of alpha diversity metrics (Chaol, ACE 
Richness, Shannon's Index, Simpson's Index) were also performed 
using MOTHUR. Richness and diversity estimates were calcu- 
lated on random subsamples set to the size of the smallest 
library to alleviate biases by sample size (Youssef and Elshahed, 
2008; Gihring et al., 2012). Percent coverage of the libraries, 
which determines the probability that all genotypes present in 
a given set of samples were recovered at least once, was calcu- 
lated as follows: [1 — (n,/N)] * 100 where n, is the number of 
unique OTU and N is the total number of clones sampled in 
the library (Good, 1953). The statistical comparison of commu- 
nity structure among the three ponds' clone libraries was tested 
using / -Libshuff and analysis of molecular variance (AM OVA) as 
implemented in Mothur software (Schloss et al, 2009). AMOVA 
tests whether the genetic diversity within communities is sig- 
nificantly different from their pooled genetic diversity (Schloss, 
2008). / -Libshuff uses the integral form of Cramer-von Mises- 
type statistic as described in Schloss et al. (2009). To account for 
multiple comparisons among the three libraries (i.e.. Ponds 9, 11, 
12), Bonferroni corrections for P- values were used to determine 
significance. 

Physico chemical data were analyzed using a Principal 
Components Analysis (PCA) using Primer software v. 6.1.11 
(Primer-E Ltd., Plymouth, UK). 

NUCLEOTIDE SEQUENCE ACCESSION NUMBERS 

Sequences were deposited in GenBank and were assigned the 
accession numbers KF234269-KF234397, KF814118-KF814651, 
and KF870833-KF870836. 



RESULTS 

PHYSICOCHEMICAL CHARACTERIZATION OF PONDS 

Total salinity increased across the three ponds with Pond 9 at 
180 g 1^' and the other two ponds being more similar at 370 
and 380 g 1^' salinity. PCA based on variation in physicochem- 
ical variables confirmed that Pond 11 and 12 were more similar 
to each other than to Pond 9 with the majority of the variation 
(98.4%) observed along the primary axis, which can be explained 
by variation in chemical species, but not temperature (Figure 1). 
Most major cations and anions increased in concentrations with 
the exceptions of Ca^+ (Table 1). This is reflected by the differ- 
ence in directionality of the eigenvector for calcium compared to 
the other ions. The drop in calcium concentration in Pond 1 1 and 
12 is likely due to precipitation of gypsum (CaS04), which begins 
in Pond 9. Somewhat surprisingly, sulfate concentrations did not 
decline in parallel with calcium concentration. The slight drop in 
sodium and chloride ions in Pond 12 water compared to Pond 1 1 
is likely due to the precipitation of halite (NaCl) at these elevated 
salinities; halite deposition was apparent in Pond 12. 

ARCHAEAL 16S rRNA SEQUENCE DIVERSITY ACROSS ESSA 
SALTERN PONDS 

Archaeal 16S rRNA sequence diversity included members of 
genera previously cultured from the ESSA saltern including 
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FIGURE 1 1 Principal Components Analysis (PCA) biplot. Vectors show 
the direction of maximum change for variables. 



Table 1 | PhysicocPiemical parameters in evaporation ponds. 
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numerous Halorubrum-\ike sequences (OTU 14-19, 90-99% 
similarity to cultured species) obtained from all three ponds at 
clonal abundances of 1-11% of the community (Figure 2). We 
also recovered one Haloarcula-like clone (98% similarity) from 
Pond 12. Pond 11 and 12 clone libraries were dominated (62 
and 70%, respectively) by a diverse assemblage of closely related 
Haloquadratum-Wke. sequences (OTU 1-6). This was especially 
true of a single lineage 99% similar to the type species H. walsbyi 
(OTU 1, 47 and 45% of Ponds 11 and 12 community, respec- 
tively). The other five, less abundant OTU in this group ranged 
from 1-12% of community and were 92-99% similar to H. wals- 
byi. Some of these lineages were highly similar (98 > 99%) to env- 
ironmental sequences obtained in Australian crystaUizer ponds 
(Oh et al., 2010) and the Santa Pola saltern in Spain (Zhaxybayeva 
et al., 2013) (Figure 2). Sequences related to Halorhabdus uta- 
hense (OTU 10) represented ~2-3% of clones in Pond 11 (98% 
similar). Additionally, lineages related to other uncultured 
haloarchaeal lineages with no known cultured representative 
within the family Halobacteriaceae were observed at all sites. 

In Pond 9, a Halorubrum-like lineage (OTU 18, 1% of 
community) and singletons (~6% total) were nearly identical 
(>99.5% similar) to environmental sequences from the Santa 
Pola saltern (Zhaxybayeva et al., 2013). Aside for these. Pond 9 
clone sequences were exclusively related to uncultured haloar- 
chaeal lineages. One lineage (OTU 12, 12% of the community) 
was >97% similar to an environmental clone from a Chinese 
saltern. However, the preponderance (79% of community) was 
a single, highly redundant lineage most closely related (99% sim- 
ilarity) to sequences recovered from the Sfax saltern in Tunisia 
(Trigui et al., 2011) (OTU 7, Figure 2). This redundancy of phy- 
lotypes in Pond 9 was reflected in much lower species richness 
and diversity metrics in this community compared with the com- 
munities in the two salt-saturated ponds (Table 2) and a flatter 
rarefaction curve in the Pond 9 sample (Figure 3A) at most sim- 
ilarity levels. We analyzed the alpha diversity results at 4 different 
OTU cut-offs: 0, 1, 3, and 5% evolutionary distance levels among 
sequences. Regardless of cut-off used, the number of OTUs cal- 
culated and the Shannon's and Simpson's diversity (shown as 
1/D = Dominance) were always lower for Pond 9 than the other 
two ponds, and the Pond 9 rarefaction curves were always lower 
than Pond 11 and 12 at each cut-off (Figure 3A). However, the 
richness estimates (Chaol and ACE) were much higher at the 0% 
cut-off than at other similarity levels within Pond 9, and both the 
estimates were even higher than Pond 1 1 at this level. These results 
suggest that there were a number of unique, but highly similar 
(<1%) OTUs in Pond 9, many in the abundant OTU 7 (data not 
shown). Significant community overlap between the archaeal 16S 
rRNA communities in Ponds 11 and 12 was confirmed with the 
/ -LIBSHUFF and AMOVA comparisons both of which showed 
no significant difference between the two libraries (Table 3), but 
highly significant differences when those libraries were compared 
with the Pond 9 library. This was also reflected in the Venn dia- 
gram showing overlap of OTUs among pond communities at the 
4 levels of similarity (Figure 3D). High degrees of overlap were 
observed between Ponds 1 1 and 12 at all similarity levels, but less 
overlap (2 OTUs) were observed between Ponds 9 and 11 and 
no overlap was observed between Ponds 9 and 12 until the 5% 



evolutionary distance threshold was employed. The overlapping 
OTUs between Ponds 9 and 1 1 were the Halorubrum-like lineage 
(OTU 18, 99% similarity) and the uncultured lineage recovered 
from a Chinese saltern (OTU12, 97% similarity). 

BACTERIORHODOPSIN GENE DIVERSITY 

We successfully created an alignment and ARB database of bop 
sequences from this study and those downloaded from the 
Genbank database and created a phylogenetic tree (Figure 4). 
Overall, similar patterns of community shifts across the ponds 
to the archaeal rRNA library were observed when bop nucleotide 
sequences were analyzed. Many environmental clones were recov- 
ered from Ponds 11 and 12 (58 and 77% of community, respec- 
tively) that had >90% bootstrap support for clustering with H. 
walsbyi on the tree (Figure 4). These varied in similarity from 
98% (OTU 4), 95% (OTUs 5-6), 90% (OTU 3) to 75% (OTUs 
1-2) to the H. walsbyi bop gene sequence. Sequences in the latter 
group were 96-100% similar to environmental clones reported 
in one of the few published environmental bop surveys from the 
Santa Pola saltern (Papke et al., 2003). Only a single sequence, 
88% similar to a cultured Halorubrum from the same study, 
was obtained from Pond 11 and no Haloarcula were found in 
the environmental bop library. This was unexpected since we 
detected both these groups with 16S rRNA genes in Pond 12 and 
we have successfully cultivated members of these genera from 
the ESSA ponds (see bolded culture bop sequences in Figure 4). 
Additionally, one Pond 11 lineage (OTU 10, Figure 4) was closely 
related (~96% similar) to environmental clones from both Santa 
Pola and a saltern in Chiku, Taiwan (Lin et al., unpublished). 
However, none of the Pond 9 clones were closely related to any 
sequences in Genbank (Figure 4). One cluster of closely related 
OTUs (7-9, 95-98% similar to each other, 71% similar to H. wals- 
byi) represented 44% of clones in Pond 9 and was only found in 
that location. Only one well-supported cluster (100% bootstrap) 
of sequences was found in both Ponds 9 and 1 1 (OTUs 12-13 plus 
singletons). This phylogroup comprised 48% of the community 
in Pond 9 and 25% in Pond 11, and again had no close cultured 
relative (70-71% similar to Natronococcus and Halobiforma bop 
sequences). As with the archaeal 16S rRNA data, higher species 
richness and diversity was found in Ponds 11 and 12 compared 
to Pond 9 (Table 2) at most OTU cutoffs, reflected in higher 
OTU redundancy in the Pond 9 rarefaction curves (Figure 3B). 
Once again the exception was the 0% evolutionary distance cut- 
off, where very high richness estimates were observed for Pond 
9 resulting in the steepest rarefaction curve of all (Figure 3B). 
For bop gene analyses, the OTU number and diversity metrics in 
Pond 9 were even higher than the other ponds at the 0% level, 
while they were the lowest for all measures with this sample at 
the 1 and 3% OTU distance level (Table 2). As with archaeal 16S 
rRNA genes, more overlap was found between Ponds 11 and 12 
at all similarity levels with only 2 OTUs (uncultured OTU 12-13) 
shared between Pond 9 and 1 1 and none between Ponds 9 and 12 
(Figure 3E). 

BACTERIAL DIVERSITY PAHERNS 

Compared with the archaeal 16S rRNA libraries, even greater 
diversity of bacterial 16S rRNA phylotypes was recovered 
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Jr Slax saltern, Tunisia clone TSLNAA9. H0157622 
UESSAPondU haloarchaeon clone 1D M13F, KF814445 
I E5SA Pond 12 haloarchaeon clone 7F, KFZ34307 OTU 1 (0/35/37) 
1 Dry Creek, Australia crystaliizer clone Cry7_14, GQ375040 
f Haloquadratum walsbyi, AY676200 
H Bajool, Australia crystallizer cione Baj_clone5, GQ374956 
ESSA Pond 12 haloarchaeon clone Z-5C_M13F, KF814410 
- Santa Pola, Spain saltern cione AAA188_I19, JN839906 
ESSA Pond 11 haloarchaeon clone Z-12E, KF814469 
Santa Pola, Spain saltern clone AAA188_E10, JNe39eB0 
. ESSA Pond 12 haloarchaeon clone 10B, KF234298 OTU 2 (0/0/10) 
ESSA Pond 11 haloarchaeon clone 3-1 □, KF234344 OTU 3 (0/5/4) 



[ v r ess; 



ESSA Pond 11 haloarchaeon clone 2-9G, KF814459 OTU 4 (0/3/G) 
ESSA Pond 11 haloarchaeon clone 2-12H, KF234334 

ESSA Pond 12 haloarchaeon clone 2'-lE M13F, KF814366 
ESSA Pond 12 haloarchaeon clone 11C_M13F, KF814261 
ESSA Pond 12 haloarchaeon clone 5B_M13F, KF814280 
Dry Creek, Australia crystallizer cione Cry7_clone 96, GQ375018 
ESSAPondll haloarchaeon clone 2-10A, KF234328 OTU 5 (0/2/0) 
ESSA Pond 11 haloarchaeon clone 3-7C partial, KF814487 
ESSA Pond 11 haloarchaeon clone 3-6C M13F, KFai4454 OTU 6 (0/1/1) 
Stax saltern, Tunisia clone 209ZH04, FN391199 

ESSA Pond 9 haloarchaeon clone 8C, KF234315 OTU 7 (83/0/0) 

Stax saltern, Tunisia cione iyHNAA17, HQ1 57576 

ESSA Pond 9 haloarchaeon clone 10C-M13F, KF814351 

- ESSA Pond 9 haloarchaeon clone 2-2H M13F, KF8143B1 

ESSA Pond 9 haloarchaeon clone 7A_M13F, KF814342 

ESSA Pond 9 haloarchaeon clone 7H M13F, KF814337 





' ESSA Pond 9 haloarchaeon clone 4C_M13F, KF814327 

Hatogeometricum borinquense. AF002984 

Sfax saltern, Tunisia cione SFGl Gill, AM947480 
ESSA Pond 11 haloarchaeon clone 3-6B, KF234318 



Dry Creek, Australia crystallizer clone Cr7_clone46, GQ374973 

- ESSA Pond 1 1 haloarchaeon clone 2-3G, KF234324 OTU 8 (0/2/0) 
Santa Pola, Spain saltern clone AAA188_F23, JNa39a90 
ESSA Pond 12 haloarchaeon clone 11D_M13F, KF814283 

ESSA Pond 11 haloarchaeon clone 3-8H, KF234327 

— ESSA Pond 11 haloarchaeon clone 3-6D, KF870833 
ESSA Pond 12 haloarchaeon clone 2-3A_M13F, KF814396 OTU 9 (0/0/2) 
Dry Creek, Auslraiia crystallizer clone Cry7_clone67, GQ375005 
ESSA Pond 12 haloarchaeon clone 10G, KF234303 

— ESSA Pond 12 haloarchaeon clone 6B, KF814306 





ESSA Pond 11 haloarchaeon clone 2-1 F, KF234319 

Sfax saltern, Tunisia clone SFHl H021 , FN391278 

ESSA Pond 11 haloarchaeon clone 2-3E, KF234323 OTU 10(0/2/0) 
Santa Pola, Spain saltern clone AAA18B_E14, JN839882 
Halorhabdus utahensis, AF071880 
ESSA Pond 12 haloarchaeon clone 7C M13F, KF614289 
Santa Pola, Spain saltern clone AAA188_N16, JN839929 
Dry Creek, Australia crystallizer done Cry7_clone50, GQ374974 

- ESSA Pond 11 haloarchaeon clone 2-1 IE, KF234330 
ESSA Pond 11 haloarchaeon clone 2-5C, KFZ34338 
ESSA Pond 9 haloarchaeon clone 2G, KF234314 OTU 11 (2/0/0) 
ESSA Pond 11 haloarchaeon clone 2-30, KF234322 OTU 12(13/2/0) 
Saifern sediment China clone YA6, EF533953 
Natronomonas pharaonis. D87971 

Haharcula japonica, EF645685 
L ^ — ESSA Pond 12 haloarchaeon clone 4B, KF234300 

■ Haharcula hispanica, AB090ie8 

J JNB39864.1, Santa Pola, Spain saltern cione AAA188_B17, JN839a64 
A- Bajool, Australia crystallizer clone Baj_cione12, GQ374959 
^ ESSA Pond 12 haloarchaeon clone 9C, KF234309 OTU 13 (0/0/2) 
Tirez lagoon, Spain cione SAlO, EU722677 
- ESSA Pond 9 haloarchaeon clone 5E, KF234313 
Santa Pola, Spain saltern clone AAA188 O07, JN839932 
— ESSA Pond 9 haloarchaeon clone 8E-M13F, KF814323 
_r-i — Dry Creek, Australia crystallizer clone Cry7_clone60, GQ375002 

ESSA Pond 12 haloarchaeon clone 4C, KF234308 OTU 14 (0/0/4) 

ESSA Pond 12 haloarchaeon clone 1G, KF234299 OTU 15 (0/2/1) 
ESSA Pond 11 haloarchaeon clone 2-11B,KF814466 OTU16(0/2/0) 



Dry Creek, Australia crystallizer clone Cry7_clone74, GO375008 
ESSA Pond 12 haloarchaeon clone 10D, KF234301 OTU 17(0/1/1) 
Sfax saltern, Tunisia clone SFG1A081, AM947495 
Halorubrum californiense, EF139B54 
ESSA Pond 12 haloarchaeon clone 2-2F_M13F, KF8i4393 



r Santa Pola, Spain saltern cione AAA188_B1S, JNS39865 

r*i ESSA Pond 9 haloarchaeon clone 10F, KF234317 0TU 18(1/8/0) 

I Halorubrum onentale. AM235786 

~ ESSA Pond 12 haloarchaeon clone 12A, KF234304 OTU 19 (0/8/2) 

Methanococcus maripaludis C5, GP000609 



■ Hsloquadratum-llke Pll/12 OTU 1 

- Hsloquadratum-llke P12 OTU 2 

- Haloquadratum-like PH/12 OTU 3 
» Haloquadratum-like Pll/12 OTU 4 

- Haloquadratum-like Pll OTU 5 
Haloquadratum-like Pll/12 OTU 6 

■Uniult P9 0TU7 
UncultPll OTU 8 

■ Halorhabdus-like PIZ OTU 9 
Halorhabdu^-like Pll OTU 10 

- Uncult P3 OTU 11 

■ Uniult P9/1I OTU 12 

■ Uncult P12 OTU H 

Ha loiu brum- like P12 OTU 14 
Halorubrum- like Pll/12 OTU 15 

■ Halorubrum-like Pll OTU 16 

■ Halorubrum-like Pll/12 OTU 17 
Halorubrum-like P09/11 OTU 18 
Halofubrum-like PI 2 OTU 19 

■ Singletons 




FIGURE 2 I Unrooted maximum likelihood dendrogram depicting 
relationships among representative archaeal 16 rRNA gene 
sequences at the 1% evolutionary difference level for clones 
obtained from the Ponds 9 (red), 11 (green), and 12 (blue) and 
closely related sequences. The tree was constructed using an alignment 
of 874 nucleotide positions (gaps and ambiguous residues were excluded 
using a custom filter in ARB). Symbols at branches represent nodes with 
bootstrap support >75% (□), 90% (o), and 100% (•) for maximum 



likelihood trees (750 replicates). Partial ESSA sequences labeled with 
M13F/R and short Santa Pola saltern sequences (accession = JN839XXX) 
were added using the parsimony addition tool in ARB. For OTU (>99% 
similarity) containing more than 1 sequence, the number of duplicates 
from each pond (9/11/12) is shown in parentheses. Methanococcus 
maripaludis was used as the outgroup. Pie charts show the relative 
representation of each of the identified OTUs for each pond, with those 
found only once (singletons) grouped. Legend is for all three charts. 
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Table 2 | 16S rRNA and bop nucleotide diversity analyses among ESSA source ponds. 



Sequence 



Archaeal 16S rRNA 



Source pond Sequences (n) OTU 



9 
11 
12 



105 
74 

83 



Coverage {%) H 



27,10,06^ 74,91,94 
3721,17 50,72,77 
42,24,16 49,71,81 



1/D 



2.0,0.8,0.6 3.6,1.6,1.4 
3.1,2.1,1.8 14.5,4.2,3.1 
3.2,2.2,1.6 171,4.6,2.5 



Ace richness Chao richness 

159,38,12 104,21,08 

103,52,26 96,29,20 

356,142,66 218,54,31 



Bop (nt) 



9 
11 

12 



66 
48 

47 



4711,05 29,84,93 
31,14,11 35,71,77 
32,16,14 32,66,70 



3.71.8,1.0 65.0,5.0,2.4 140,33,06 10714,05 
3.2,2.2,2.0 31.3,8.2,6.6 134,28,24 64,19,21 
3.3,2.0,1.8 30.0,5.1,4.0 15766,80 132,49,42 



Bacterial 16S rRNA 9 
11 
12 



81 
69 
77 



39,18,16 52,78,80 
170707 75,90,90 
4731,24 39,60,69 



3.2.2.2.2.0 175,6.3,5.5 222,50,43 141,30,25 

1.8.1.1.1.1 3.4,2.0,2.0 45,0707 32,0707 
3.6,3.0,2.5 39.6,16.6,9.1 386,184,370 399,108,69 



H, Shannon-Weaver Diversity Index; D, Simpson's diversity index; shown as 7/D = Dominance. 

^ Commas separate mean estimates calculated using 0, ?, and 3% OTU dissimilarity cut-offs respectively for all analyses. 



from the ESSA ponds. This included members of the Alpha-, 
Delta-, and Gammaproteobacteria, Bacteroidetes, Firmicutes, 
Verrucomicrobia, and algal plastid sequences (Figures). 
Pond 9 communities had relatively high representation of 
Verrucomicrobial Puniceicoccus-\ike sequences (OTU 10-11, 
22% of clone), but the largest proportion (47%) was found in 
two Gammaproteobacterial phylogroups (Figures). The first 
of these (OTU 12, 14% of community) was 97% similar to 
Spiribacter salinus, a recently isolated photoheterotroph that was 
found to be abundant (~16% of community) at ~19% salinity 
in the Santa Pola saltern in Spain (Ghai et al., 2011; Leon et al, 
2013). The second (33% of community) was 99% similar to an 
environmental clone from the Salton Sea, a modestly hypersaline 
(~40 g 1^^) lake in southern California. 

Pond 11 diversity was quite limited, not only in diver- 
sity metrics, but also taxonomically since only members of 
the Phylum Bacteroidetes were detected, including relatives of 
Psychroflexus, Sediminibacterium, Owenweeksia, and Salinibacter 
(Figures). Psychroflexus-like (~90% similarity to cultured 
species) sequences comprised ~16% of the Pond 11 library, 
Sediminibacterium-\ike species (94% similarity to S. salmoneum) 
made up over 7% of the community. Only three sequences (2 
from Pond 11,1 from Pond 12) were closely related (99% similar) 
to a Salinibacter culture isolated from ESSA (Sabet et al., 2009), 
which is >98% similar to the type species S. ruber. However, 
69% of the Pond 1 1 library was comprised of a single uncultured 
Bacteroidetes lineage (OTU 5, <80% similar to nearest cultured 
representative Saprospira grandis) . This lineage was closely related 
(>98% similar) to environmental clones recovered from Lake 
Tuz in Turkey and a saltern in Chula Vista, CA (Zhaxybayeva 
et al., 2013). This OTU also made up ~14% of the community 
in Pond 12. Despite overlap in this phylogroup. Pond 12 differed 
from Pond 1 1 in having a greater number of Deltaproteobacterial 
sequences (40% of community). This included one abundant 
group (13% of community) that was <85% similar to cultured 
sulfate-reducing species such as Desulfobacca acetoxidans, but 
was 96-99% similar to metagenomic sequences recently obtained 
in Lake TyreU in Australia (Podell et al., 2013). Chlorophyte 
plastid sequences comprised 10% of the Pond 12 community. 



This included one sequence related to the halophUic green alga, 
Dunaliella (~97% similar), but most were only distantly related 
to cultured algae (<80% similar to Monomastix sequences). 
Nearly one third of sequences in Pond 12 were singletons, found 
only once in the library. 

The low number of phylotypes in Pond 11 was reflected in 
low richness and diversity estimates (Table 2) and a nearly flat 
rarefaction curves for the Pond 11 community (Figure 3C). By 
contrast, the diversity and richness values in Ponds 9 and 12 were 
much higher. In Pond 9, there was less OTU redundancy for bac- 
teria than the archaeal community. Somewhat surprisingly, at all 
OTU levels, richness and diversity metrics were highest in Pond 
12 compared to the other ponds, exceeding Shannon index val- 
ues of 3.0 at the 0% evolutionary distance. In contrast to the 
archaeal 16S rRNA results, there was almost no overlap among 
bacterial sequences obtained from the three ponds with only 
2 shared OTUs (OTU 1,5) between communities (Figure 3F). 
This was confirmed by pairwise / -LIBSHUFF and AMOVA com- 
parisons that revealed highly significant differences between all 
ponds including 11 and 12 (Table 3). 

DISCUSSION 
ARCHAEAL DIVERSITY 

Archaeal 16S rRNA sequence diversity estimates in Ponds 11 and 
12 were comparable to those reported for the Santa Pola saltern 
and somewhat lower than in the Sfax ponds (Baati et al, 2008). 
The increase in archaeal diversity metrics at higher salinities was 
driven primarily by increased representation of culturable haloar- 
chaeal groups including Halorubrum and Haloarcula, groups that 
have been previously isolated from the ESSA saltern (Sabet et al., 
2009). Over 60% of sequences recovered in these two ponds were 
Haloquadratum-like. Among these, the most abundant lineage 
recovered from Ponds 11 and 12 (>40% in each) was only 1% 
divergent from the type species Haloquadratum walsbyi. We also 
identified sequences with lower similarity (93-97%) to H. wals- 
byi, that clustered in the phylogeny with a clone recovered in a 
more recent study in Australian crystallizers (Oh et al., 2010). 
We suggest that these are likely members of another species of 
Haloquadratum that has yet to be cultivated. Variation among 
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Sequences Sampled 




Sequences Sampled 

P9 - 0% P9-1% — — P9-3% P9-5'/. 

P11 - 0% ^— P11-1% — — P11-3% P11-5% 

P12-0% P12-1% — — P12-3% P12-5% 

FIGURE 3 I Rarefaction curves and Venn diagrams showing the shared in red, Pond 11 in green, and Pond 12 in blue. Analyses were performed using 
and unique OTUs across Ponds for archaeal 16S rRNA genes (A,D), bop a 0, 1, 3, 5% evolutionary distance cut-offs shown in legend for rarefaction 
genes (B,E), and bacterial 16S rRNA genes (C,F). Pond 9 results are shown and in parentheses (0,1,3,5) for Venn. 



Haloquadratum-\ike genotypes assayed by DGGE in different 
ponds in the Sfax saltern has also been reported (Boujelben 
et al., 2012). It has been speculated that the seemingly cos- 
mopolitan occurrence of Haloquadratum in the highest salinity 
ponds in salterns may be due to their tolerance of extreme salin- 
ity fluctuations and dispersal mechanisms (Oh et al., 2010), 
although at this point little data is available especially regarding 
dispersal. 



The haloarchaeal community was much less diverse in Pond 
9, the main site of gypsum precipitation (CaS04) in the ESSA 
saltern. The Pond 9 archaeal library was dominated by uncultured 
lineages, especially a single phylotype >99% identical to an envi- 
ronmental sequence recovered from the Sfax multipond saltern 
in Tunisia (Baati et al., 2010). Interestingly, the phylotype from 
that study was also recovered from a pond with similar salinity 
(~18%), suggesting this group maybe selected by this salinity or 
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Table 3 | Community comparisons among ponds for 16S rRNA 
sequence libraries. 



Library 


Statistic 


Comparison 


P-value 


Archaea 


AMOVA'i 


9 vs. 11 
9 vs. 12 

11 vs. 12 


<0.001 
<0.001 
0.22 




/-Libshuffb 


9 vs. 11 
9 vs. 12 

11 vs. 12 


<0.001 (XY), <0.001 (YX) 
<0.001 (XY), <0.001 (YX) 
0.0002 (XY), 0.1776 (YX) 


Bacteria 


AM OVA 


9 vs. 11 
9 vs. 12 
11 vs. 12 


<0.001 
<0.001 
<0.001 




/-Libshuff 


9 vs. 11 
9 vs. 12 
11 vs. 12 


<0.001 (XY), <0.001 (YX) 
<0.001 (XY), <0.001 (YX) 
<0.001 (XY), <0.001 (YX) 



''ForAMOVA analyses, P < 0.017 was considered significant based on Bonferroni 
corrections for multiple comparisons (significant comparisons bolded). 
^For Libshuff analyses, pairwise comparisons were made in both directions, so 
P < 0. 0085 in both directions was considered significant (in bold) based on the 
Bonferroni correction (Singleton et al., 2001). 



some other physicochemical factors across these geographically 
isolated coastal salterns. 

Bop GENE DIVERSITY 

Our findings from archaeal ribosomal gene sequences were largely 
confirmed in our analyses of the functional gene (bop) coding 
for bacteriorhodopsins. Overall, we found a diverse assemblage 
of sequences that clustered within the Halobacterales with com- 
parable rhodopsin gene diversity (H' = 1.0-2.2 at the 1 and 
3% OTU level) compared with past studies of bop gene diver- 
sity (H' = ~1.4) (Papke et al, 2003; Pasic et al, 2005). There 
was greater overlap in the sequences recovered in ESSA with 
those from the Santa Pola saltern (Papke et al, 2003) than the 
Seeovlje saltern, which was not found to have Haloquadratutn- 
like sequences (Pasic et al., 2005). However, in contrast with Papke 
et al. (2003), we found more evidence of potential pond-specific 
communities (i.e., bop phylotypes only found in one of the three 
ponds). A single bop gene lineage was found to have similar- 
ity with a database sequence from a saltern in Chiku, Taiwan 
(Lin et al, unpublished). No overlap was observed with bacte- 
riorhodopsin sequences recently reported from the Dead sea, an 
athalassohaline habitat (Bodaker et al., 2012). 

Many of the sequences we recovered, especially in Pond 9, were 
from uncultured lineages, neither closely related to any in the 
NCBI database, nor closely related to bop sequences from cul- 
tures recovered in the ESSA saltern. The sequence novelty may be 
due in part to unique ESSA-specific phylotypes, but is also likely 
due to global undersampling and the relative paucity of available 
environmental bop gene sequence data in databases. 

Since this is one of the only studies with directly parallel 
16S rRNA and bop gene sequencing from the environment, we 
can attempt to overcome some of the limitations of the limited 
database for identification. For example, in the archaeal 16S rRNA 



library in Pond 9, a single phylotype was highly abundant (79% 
of library) and exclusively found in this pond. In our bop gene 
library, an abundant cluster of closely related OTUs (44% of 
library) was also exclusively found in Pond 9. These sequences 
may all represent variants of bacteriorhodopsin within the same 
species, as functional genes typically display more genetic varia- 
tion than highly conserved 16S rRNA gene. 

BACTERIAL DIVERSITY 

Among our ESSA Pond 9 bacterial sequences, the largest group 
recovered was 99% similar to an uncultured gammaproteobac- 
terial sequence from the Salton Sea, a moderately saline (~40g 
1^^) endorheic lake in southern California (Dillon et al., 2009a). 
This suggests that this lineage may live within the lower range of 
hypersaline conditions and may explain why it was not recovered 
in Ponds 1 1 and 12. Interestingly, a number of bacterial sequences 
from Pond 9 were closely related (~95% similar) to sequences 
recovered in a previous study of the photosynthetic microbial 
mats found in Pond 4 of the ESSA saltern (Ley et al., 2006) as well 
as evaporitic mats from EUat, Israel (>99% similarity) (Sorensen 
et al, 2005). This suggests that in interconnected saltern sys- 
tems, the microbial mats found at lower salinities may serve as 
a source of bacteria resident in the plankton further up the salin- 
ity gradient, although recovery of genes via PCR-based methods 
does not confirm that they were actively growing at these higher 
salinities. 

Algal plastid sequences were obtained in Pond 12. We recov- 
ered a single clone of the halophilic alga Dunaliella, which was 
somewhat unexpected since a past report noted the absence of this 
group from the ESSA ponds despite its prevalence in other more 
nutrient-rich salterns (lavor, 1983). However, the most abun- 
dant group of plastids sequences (10% of clones in Pond 12) 
showed only modest (~94%) similarity with an environmental 
clone found in Lake Tebenquiche in Chile (Demergasso et al., 
2008), and among cultured relatives was distantly related (<80% 
similar) to the freshwater Prymnesiophyte Monomastix (Turmel 
et al., 2009), suggesting this may represent a previously unknown 
halophilic algal lineage. 

In contrast to the other two ponds, which had relatively 
high numbers of taxonomic groups, only members of the 
Phylum Bacteroidetes were recovered in Pond 11. We found high 
clonal abundance (~27% of community) of sequences related to 
Psychroflexus, Owenweeksia, and Sediminibacterium. Psychroflexus 
strains have been previously cultured from hypersaline habitats 
(Donachie et al., 2004; Zhang et al., 2010) and members of this 
genus have been observed in salterns and high altitude athalasso- 
haline lakes in Tibet and Chile (Benlloch et al., 2002; Wu et al., 
2006; Dorador et al., 2009). No known halophilic members have 
been cultured from the Owenweeksia and Sediminibacterium gen- 
era, although members of these lineages have been isolated from 
marine habitats (Lau et al, 2005; Khan et al., 2007). These find- 
ings contrast with past studies of hypersaline lakes and salterns 
where members of the Salinibacter genus were abundant (Anton 
et al., 2000; Benlloch et al, 2002; Demergasso et al, 2004; Baati 
et al., 2008). We only rarely detected Salinibacter-like lineages in 
ESSA, with two clones from Pond 11 closely related (>99% sim- 
ilar) to a cultured ESSA Salinibacter (Sabet et al, 2009) and one 
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I Santa Pola, Spain bop clone, AY366560 

I ESSA Pond 12 bop clone D8, KF234280 OTU 1 (0/3/15) 
r^ESSA Pond 12 bop clone C2, KF814131 OTU 2 (0/2/3) 
ESSA Pond 11 bop clone BE, KF234296 



rV 



ESSA Pond 11 bop clone 6H, KF234295 OTU 3 (0/10/15) 
ESSA Pond 12 bop clone E7, KF234282 
■ ESSA Pond 11 bop clone 9C, KF234297 OTU 4 (0/1 0/1) 

Haloquadratum walsbyi, AM 180088 

ESSA Pond 11 bop clone 11A, KF814217 0TU 5 (0/3/0) 
ESSA Pond 12 bop clone C6, ARB_3B2C9B8E OTU 6 (0/0/2) 
Santa Pola, Spain bop clone, AY366555 
ESSA Pond 11 bop clone 6A, KF234293 
ESSA Pond 12 bop clone B8, KF234276 
- ESSA Pond 12 bop clone E12, KF234281 



▼ E5 



ESSA Pond 12 bop clone A2, KF234271 



ESSA Pond 9 bop clone 2-D2, KF234287 OTU 7 (9/0/0) 
ESSA Pond 9 bop clone 2-D3, KF814194 

ESSA Pond 9 bop clone 2-A9, KF814174 OTU 8 (2/0/0) 
ESSA Pond 9 bop clone 2-B12, KF814176 
ESSA Pond 9 bop clone 2-C9, KF81 41 89 OTU 9 (1 8/0/0) 
Santa Pola, Spain Fosmid FLAS10H9, EF558547 
ESSA Pond 11 bop clone 6C, KF234294 
ESSA Pond 9 bop clone C8, KF234284 

^ cultured ESSA Halorubrum isolate RBC-B-10-2, KF814254 
JT Santa Pola, Spain Halorubrum sp., AM777510 
pT ^ cultured ESSA Halorubrum isolate 11GM-10-3, KF814255 

i-| ^ Santa Pola, Spain Halorubrum sp., EF558558 

Santa Pola, Spain Halorubrum sp., AM777540 

Santa Pola, Spain Halorubrum sp-, AM777438 
Santa Pola, Spain Halorubrum sp., AM7774S5 
Santa Pola, Spam Halorubrum sp., A(V1777433 
Santa Pola, Spain Halorubrum sp,, AM777427 

ESSA Pond 11 bop clone 3A, KF234292 
/*. Chiku, Taiwan clone, HM4751 1 3 
Santa Pola, Spain clone, AY366553 
~^ ESSA Pond 11 bop clone 118, KF234289 OTU10(0/3/0) 
Chiku, Taiwan clone. HM475120 




4^ 



ESSA Pond 12 bop clone A7,KF234272 

cultured ESSA Haloarcula isolate 9-D-U, KF814252 
unknown location Haloarcula sp., AY279550 
Yiwu Lake, China Haloarcula sp., D0146999 
T ^ cultured ESSA Haloarcula Isolate RBC-A-10-2, KF814250 
^ cultured ESSA Haloarcula isolate 9-B-U, KF814251 

Chiku. Taiwan clone. HM475129 

cultured ESSA isolate Bitterns-U, KF814253 
Santa Pola, Spain isolate TP072, EF558555 
Santa Pola, Spain isolate TP100, EF558556 
— ESSA Pond 11 bop clone 2B, KF234291 

Halosimplex carlsbadense. EF558554 
ESSA Pond 12 bop clone B2, KF234273 



ESSA Pond 12 bop clone C11, KF234277 

ESSA Pond 9 bop clone D4, KF234285 OTU11 (2/0/0) 

J ESSA Pond 12 bop clone B3, KF234274 

1- ESSA Pond 12 bop clone A11, KF234269 
ESSA Pond 12 bop clone C9, KF234278 
ESSA Pond 12 bop clone D1, KF234279 
ESSA Pond 9 bop plate 2 clone E9, KF814210 OTU 12 (10/2/0) 
ESSA Pond 11 bop clone 3C, KF814227 
ESSA Pond 9 bop clone F5, KF814165 
ESSA Pond 9 bop clone 2-F2, KF814212 
ESSA Pond 9 bop clone C5, KF234283 OTU 13 (20/9/0) 



Halobiforma /ac/sa/si, AY279551 
- Natronococcus aibiensis, AY27954 



P9 



P11 



P12 



■ Uncult Pll/12 OTU 1 

■ UncultPll/12 0TU Z 

- Haloquadratum-like Pll/12 OTU 3 

■ Hainquadratum-like Pll/lZ OTU 4 

- Haloquadratum-like Pll OTU 5 
Haloquadratum-lik^ P12 OTU 6 
Uncult P9 OTU 7 

■ Uncult P9 OTU 8 
Uncult P9 OTU 9 

■ Uncult Pll OTU 10 
> Uncult P9 OTU 11 

Uncult P9/11 OTU 12 
Uncult P9/11 OTU 13 

■ Singletons 




FIGURE 4 I Unrooted maximum likelihood dendrogram depicting 
relationships among representative bop gene sequences at the 1% 
evolutionary difference level for clones obtained from the Ponds 9 (red), 
11 (green), and 12 (blue), ESSA cultures (bold) and closely related 
sequences. Tlie tree was constructed using an alignment of 312 nucleotide 
positions (gaps and ambiguous residues were excluded using a custom filter 
in ARB). Symbols at branches represent nodes with bootstrap support >75% 



(□), 90% (o), and 100% (•) for maximum likelihood trees (550 replicates). For 
OTU (>99% similarity) containing more than 1 sequence, the number of 
duplicates from each pond (9/11/12) is shown in parentheses. Halobiforma 
lacisaisi and Natronococcus aibiensis sequences were used as outgroups. 
Pie charts show the relative representation of each of the identified OTUs for 
each pond, with those found only once (singletons) grouped. Legend is for all 
three charts. 
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Lake TsbarquicHe, CHila Clors ATlBOctCS, AVB62745 
E5SA Pond 1 1 bacterium cJone 3E, KF234390 OTU 1 (TJGIO) 
ESSA, GJerrern Negm, MX Evaporile Clone 1Y25, EF1Q59eB 

Pond 11 bacterium clone 5fl, KFZ34386 OTU 2 (OIM) 
Stax saltern, Tunisia dore, CU467«0 

Psych roflBKus Ep,. EF 190063 
PsychfoHe/uB sEdmiins. EU135715 
PsydifoneKUS tiopicus. AF51 3434 
ESSA Pond 9 bacterium clone 3C, KF234351 
endoevaporllic mal clone E4aFl I, DQ103639 
ESSA Pond 9 bacterium clone 2A. KF23435T OTU 3 (7/0/0) 
ESSA Pond 11 bacterium clone 7A, KF234385 0TU4 (0/2/0) 
Owein/teeksia hongkongensis, ABIESOea 

■.lone HS-EN62B, FJ536453 
ESSA Pond 12 bacterium clone 11H, KFZ34361 
irael endoavapomic mat Clone E6aC02, DO103645 
Israel endQevaporitlc mat clone E6aEl1. DO103643 
ESSA Pond 12 bacterium clone 11B, KFZ343eO 
Cabn Rojo, Pueno Rko Clone MAT-CR-H6-DQ2. EU245346 

" " SA sallBcn clone AAA186_B06, JNB39e54 
.Ione4-01B. EF45S715 

actcrium clone 3H, KFZ34392 OTU 5 [o;48/1 1) 
aclerium clone 7G_M13R, KFB14633 

Chile done A1120ctB3, A¥a62793 
ESSA Pond 12 bacterium clone 1D, KF2343T4 OTU6(0/Q/5) 
endoevaponlic mal clone E28B05. D01Q364B 
baled reactor clone 35, FJ623301 

ESSA Pond 1 1 bacterium clone 6E, KF2343gi OTU 7 (0/5/0) 

Sediminibaulenum salmoneurr, EF407879 
Flai/isolibaclsfginsengiSeiree, AB267476 



Pond 12 bacterium clone 10F, KF2343E9 

PI GOS malne melagenome, AACY020E25321 

JF Salton Ssa. CA. USA clone SSI/V6IAp, E 

^ ESSA Pond 9 bacterium clone GB, KF234354 OTU S (2/0/0} 
EF9aB655 

ESSA Pond 9 bacterium clone 2-60-M13F, KFS14S79 
bacterium clone 10O.ft113F, KFS14531 
ESSA Pond 12 bacterium clone BH, KF234370 

EU24B149 

ESSA Pond 12 bacterium clone SE, KF2343eS 

ap. CH-IO'S.EUaSIBII 
ESSA Pond 1 1 bacterium clone 4C , KFZ34384 OTU 9 (0/2/0) 
I, FN39431 
uncultured Saffn/bacler sp . AYe62783 
8_D16 

ESSA Pond 12 bacterium clone BD, 

Ameur Lake, Algeria clone ARDBACSS2, EU8S9375 
Salmibsaer iiiber. CP000159 
Pijiiiceicoccus i/aimicol3, DO539046 

GN01-0.0004, 00154337 
ESSA Pond 9 bacterium clone ZF M13F, KFS14S14 OTU 10 (13/0/0) 
ESSA Pond 9 bacterium clone 3E, KF23439S OTU 11 (5/0/0) 
Pond 9 bacterium clone Z-SG.MtSF, KF814592 
ESSA Pond 12 bacterium clone 2E_M1 3F, KF870B34 
ESSA Pond 12 bacterium clone 5C, KF2343S3 

Its clone HS-EN2B, FJ536428 

VKlivallts i/adensis. AY049713 
Sfax sallem, Tunisia clone 3FDia071, FN393476 

Pond 12 bacterium clone 1ZG, KFZ34362 
Halorhodospira hahpMa. CP000544 
Jidong OiHI Bid. China clone SN1 13. Eli735eBl 

ESSA Pond 12 bi 

niofiaiospira haiophtla. □Q469576 

NR_103B52 

CU 4674 60 

ESSA Pond 9 bacterium clone 8G, KF234397 OTU 12 (11/010) 
Salton Sea. CA, USA clone SSWaSAu, EU592377 

9 bacterium clone TB, KF234396 OTU 13(27/0/0) 

Tunisia clone, CU467463 

AJ64D132 

ESSA Pond 9 bacterium clone 2C, KFZ3435S 

Saliniv'biio bud3:i. EF583S83 
Aeromonas hydiophila. CP0004ea 
marina W-JKIO. EF554894 
;e Teicoco, MX soil dons TX4CB_93, FJ152962 
Pond 12 bacterium clone 5D, KF234364 
RickBltsiellagiylli, U97547 

ESSA Pond 12 bacterium clone 5H, KFai4547 OTU 14 (0/0/2) 
eNil5Aij0613_710, FJ638317 

deniliificans. D0469582 
Rop, Puerto Rico mat clone MAT-CR-P2-C0B. eU24604B 
Pond 12 bacterium clone 4G, KFZ343BZ 

;ala' de Ascotan. Chile done Asc-3-74, EF632563 
ESSA Pond 9 bacterium clone 3G, KF2343B2 



RhodoBacteraceae bacterium ROf 
ESSA Pond 9 bacterium cloi 

Roseovarius pacrticus, DO120726 



, Aya41782 

! 9D M13F, KF814511 



ESSA Pond 9 bacterium clone 1G, KF234356 CTU 15(3/0/0) 
ESSA Pond 9 bacterium clone 3A M13F, KFB14497 

Linflsey Lake. NM, USA dona LL39b, EF1056B9 
ESSA Pond 12 bacterium clone SF, KF234365 

Rhodovibdo saf/n; 

Sidi Ameur Lake, Algeria clone ARDBACSS20, EUE693a2 
Pond 12 bacterium clone 4A, KF234379 
Lake. Algeria ctone ARDBACSS44. EU8S93e6 
ESSA Pond 12 bacterium clone 4B, KF234380 
ESSA Pond 12 bacterium clone 1C, KFZ34373 
ESSA Pond 12 bacterium clone 10B M13F, KFB14529 

(O/Q/5) 

(0/0/10) 

lac78_Flocs. AB49ia4G 
SOU dona p27knok, FJ4 79048 
Pond 9 bacterium clone 5C, KF234353 

EU220836 

ESSA Pond 12 bacterium clone 4 E, KF8i4S55 OTU18(0/Q/6| 

12 bacterium clone 8A, KF2343B3 
ESSA Pond 12 bacterium clone 3F,KF814554 OTU 19 (0/0/7) 
Stajc saltetn, Tunisia clone SFD1G071. FN3934aa 
DEsdfohalabiiim u/a/iense. DQ067421 




clone 1A, KF2343T2 
I clone 3B, KF234377 OTU 20 (0/0/8) 
:lone 4D, KFa7Qa36 



Desullovibno hslophik 

ESSA 

ESSA Pond 12baeteri 
ESSA Pond 12 bacterial 

sp. OKE-1, FJ4934B7 
^ Lake Tebanqulche, Chile Dunaliella clone AliaAugA12, 
' — ESSA Pond 9 bacterium clone 2-6E M13F, KF370B35 

Synec/rocys/iH Bp FCC B702, ABD4153B 

^ Cabo Rojo, Puerto Rco mal done MAT-CB-H5-E08. EU245295 
^ ESSA Pond 12 bacterium clone 9E, KF2343T1 

— Halaneerobium salsuginis. L22B90 

ESSA Pond 9 bacterium clone 2-4C M13F, KFB14564 OTU 21 (3/0/0) 
FB04D02, FM672896 
ESSA Pond 12 bacterium dona 3C, KF234393 

VF5. AE000657 



FIGURE 5 I Unrooted maximum likelihood dendrogram depicting 
relationships among representative bacterial 16 rRNA gene 
sequences at the 1% evolutionary difference level for clones 
obtained from the Ponds 9 (red), 11 (green), and 12 (blue) and 
closely related sequences. The tree was constructed using an alignment 
of 1033 nucleotide positions (gaps and ambiguous residues were 
excluded using a custom filter in ARB). Symbols at branches represent 
nodes with bootstrap support >75% {□), 90% (o), and 100% (•) for 



maximum likelihood trees (250 replicates). Partial sequences labeled with 
M13F/R and short Chula Vista saltern sequences (accession = 
JN839XXX) were added using the parsimony addition tool in ARB. For 
OTU (>99% similarity) containing more than 1 sequence, the number of 
duplicates from each pond (9/11/12) is shown in parentheses. Aquifex 
aeolicus was used as the outgroup. Pie charts show the relative 
representation of each of the identified OTUs for each pond, with those 
found only once (singletons) grouped. Legend is for all three charts. 
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clone 98% similar to an environmental clone from the saltern in 
Chula Vista, CA, USA (Zhaxybayeva et al, 2013). 

Additionally, an uncultured Bacteroidetes lineage >98% sim- 
ilar to clones from Lake Tuz in Turkey (Mutlu et al., 2008) and 
the Chula Vista, CA, USA saltern (Zhaxybayeva et al., 2013) com- 
prised 69 and 14% of all bacterial 16S rRNA clones recovered 
from Ponds 11 and 12, respectively. This group was also the most 
commonly identified (29/58 bacteria) using single cell genomics 
approaches from high salinity ponds (320-350 g 1^^ ) in the Chula 
Vista saltern, which is near the US-Mexico border. Both the Chula 
Vista and ESSA salterns are derived from evaporation of pacific 
coastal waters and are less than 500 miles apart. The Zhaxybayeva 
et al. (2013) study found little overlap in bacterial communi- 
ties between the Chula Vista saltern and the Santa Pola saltern 
in Spain, suggesting geographic or other unidentified environ- 
mental differences may be responsible. The abundance of this 
Bacteroidetes phylotype in Pond 11 of ESSA and Chula Vista 
at high salinity suggests there may be environmental factors in 
common or dispersal mechanisms between the two salterns that 
explains this. However, not all groups showed the same pattern. 
Gammaproteobacteria from commonly cultivated genera such as 
Salicola and Halomonas sequences were abundant in the Chula 
Vista saltern. Members of these bacterial genera were not detected 
in this study, despite being isolated from Ponds 9 and 1 1 in a culti- 
vation study performed using samples collected in parallel to this 
one (Sabet et al, 2009). The difference in abundance of these cul- 
turable groups between the ESSA and Chula Vista salterns may be 
due to differences in nutrient levels as the latter saltern has been 
reported to be eutrophic (Javor, 1989). Elevated nutrients in the 
Chula Vista ponds may more closely resemble culture conditions 
that favor those lineages. 

The novel bacterial community members in Pond 12 were 
primarily Deltaproteobacteria. Deltaproteobacteria, especially 
sulfate-reducing lineages, have been commonly identified in 
saltern sediments (Baati et al., 2010; Lopez-Lopez et al, 2010) and 
benthic mats (Caumette et al., 1994; Four(;ans et al., 2004; Dillon 
et al, 2009b), but less commonly reported in saltern waters. We 
would not expect the lineages we recovered in these well-mixed, 
aerobic ponds to be sulfate reducers. One abundant cluster (17% 
of community) was not closely related to cultured members of 
this subphylum (<15% similar to cultured Deltaproteobacterial 
species), but was 96-99% similar to metagenomic sequences 
recently obtained in Lake Tyrell in Australia (Podell et al., 2013) 
suggesting it may be globally distributed. Overall, these findings 
combined with other recent studies (liang et al., 2007; Pagaling 
et al., 2009; Baati et al., 2010) indicate that Salinibacter may not 
always be the most abundant bacterial type in saturated brines 
and that other Bacteroidetes as well as Proteobacteria may be 
similarly well-adapted to such extreme salinities and should be 
targeted for cultivation and further study. 

DIVERSITY ALONG THE SALINITY GRADIENT 

Our findings seemed to follow our hypothesized pattern of 
increasing archaeal and declining bacterial diversity along the 
salinity gradient from Ponds 9 and 11, but the increase in bac- 
terial diversity, even higher than archaeal diversity, in Pond 12 
runs counter to this. This was surprising, since most studies in 



salterns have found that bacterial diversity declines with increas- 
ing salinity. For example, studies comparing the relative abun- 
dance of bacterial and archaeal 16S rRNA clones in the Santa 
Pola and Sfax multipond salterns found very low bacterial diver- 
sity (i.e., H' < 1.0) above 300 g salinity (Benlloch et al, 2002; 
Casamayor et al., 2002; Baati et al., 2008). Bacterial phylotypes 
have been found to outnumber archaea in an athalassohaline 
lake in the Atacama desert (Demergasso et al, 2004) and an 
alkaline, hypersaline depression in the Sahara (Mesbah et al., 
2007), but the dramatic differences in bacterial populations we 
observed between Ponds 11 and 12 were unexpected given the 
similar chemical nature of these two ponds. The only obvious 
difference between the two sites was the higher degree of precipi- 
tation of halite in Pond 12, although this does not rule out some 
unmeasured physicochemical difference between the two ponds. 

Of course, it must be noted that what we were measuring in 
this study (like the majority of similar studies) is the relative abun- 
dance of phylotypes recovered using PCR on a limited sample set, 
not direct environmental abundances. We used different primer 
sets for the bacterial and archaeal clone libraries and did not 
use quantitative PCR, so we cannot directly compare the clonal 
abundance of bacteria and archaea in these ponds and we cannot 
assume that the relative clonal abundance represents the actual 
abundance of cells in the ponds. These patterns are intriguing, but 
must be confirmed with more quantitative methods (e.g., FISH). 

What is clear is that even though the evaporation ponds in 
the ESSA multipond salterns are interconnected, with evaporat- 
ing seawater being pumped between the ponds, there are unique, 
diverse communities of both bacteria and archaea, including 
diverse bacteriorhodopsin-containing lineages, found in each. 
Future studies using metagenomic approaches to connect func- 
tional genes with taxonomic genes and targeted cultivation of 
abundant lineages identified in this study are warranted. 
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